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ABSTRACT 

We use a large, heterogeneous sample of local active galactic nuclei (AGN) that includes Seyfert Is, 
Seyfert 2s and PG quasars to investigate for the first time the relation between black hole mass (Mbh) 
and mid-infrared nuclear emission. We find a clear relation between Mbh and 10 /xm nuclear luminosity 
for these local AGNs. There are no significant differences between type 1 and type 2 objects, implying 
that the reprocessing of the 10 /tin nuclear emission is not severely affected by geometric and optical 
depth effects. We also confirm that Mbh is related to the 2 — 10 keV X-ray luminosity, but only for the 
Compton thin galaxies. We present a theoretical basis for these empirical relations and discuss possible 
reasons for the observed scatter. Our results show that rest-frame 10 /xm and hard X-ray luminosities 
(especially the former, which is applicable to all AGN types) can be powerful tools for conducting a 
census of BH masses at high redshift and for probing their cosmological evolution. 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: nuclei - quasars: general - black hole 
physics 



1. INTRODUCTION 

The cosmological evolution of the black hole (BH) pop- 
ulation is a fundamental property to constrain models of 
galaxy formation and evolution. In the last few years we 
have seen significant progress, especially for local galaxies 
where BH masses are now being measured. More impor- 
tantly, the discovery of two fundamental correlations be- 
tween the BH mass and the bulge luminosity (Kormendy 
& Richstone 1995; Magorrian et al. 1998) and the BH 
mass and the stellar velocity dispersion (Gebhardt et al. 
2000; Ferrarese & Merritt 2000) suggests a tight link in 
the evolution of the BH and its galactic host. 

To produce reliable statistics of BH masses at high red- 
shift, one would need unrealistic amounts of observing 
time if the methods employed to determine BH masses 
in the Local Universe (see Section 2) were to be used. 
Thus, understanding the cosmological evolution of BHs re- 
quires finding a good indicator of the active galactic nuclei 
(AGN) luminosity, independent from the effects of obscu- 
ration and the viewing angle, and directly related to the 
BH mass. 

Among the suggested AGN power indicators are the in- 
frared luminosity, the narrow line region [Oiii]A5007 line 
emission, the hard X-ray (2 — 10 keV) emission, and the ra- 
dio emission (e.g., Mulchaey et al. 1994). The 2 - lOkeV 
hard X-ray emission is only a good indicator of the in- 
trinsic luminosity of the AGN for those cases where it is 
transmitted through the torus, that is, in Compton thin 
galaxies. 

QSOs emit a significant fraction of their bolometric lu- 
minosity in the infrared (e.g., Sanders et al. 1989). The 
mid- and far-infrared emission in Sys and radio quiet 
quasars is predominantly thermal in origin, namely dust 
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emission (e.g., Rieke 1978; Barvainis 1987; McAlary & 
Rieke 1988; Sanders et al. 1989; and recently Haas et 
al. 2000; Polletta et al. 2000). Moreover, Spinoglio, & 
Malkan (1989) found that the 12 /xm flux is approximately 
a constant fraction of the bolometric flux in Seyfert (Sy) 
galaxies and quasars. These considerations suggest that 
the main parameter governing the mid-infrared behavior 
is the power of the central engine, almost independently 
of the geometry, dust content and dust properties of the 
central region. 

For type 2 objects, the infrared emission at A <J 5 /xm 
may still be significantly affected by viewing angle ef- 
fects (obscuration) and/or stellar contribution (Alonso- 
Herrero et al. 2001 and references therein). At A ;> 20 /xm 
there may be significant contributions from star forma- 
tion and/or the underlying galaxy. These contributions 
are not easily disentangled from the AGN emission unless 
high spatial resolution imaging is used. Thus the nuclear 
10 /xm luminosity appears as a good choice to represent 
the AGN bolometric luminosity. 

In this letter we explore for the first time the relation 
between the BH mass and the thermal 10 /xm nuclear emis- 
sion for a heterogeneous sample of low-redshift AGN: Sy 
Is, Sy 2s, and PG quasars, and explore the possibility of 
using the mid-infrared luminosity of AGN as a represen- 
tation of their BH mass. We also discuss implications for 
the BH growth and demographics over cosmological times. 

2. THE DATA AND CORRELATIONS 

We have compiled BH masses of local AGNs - Sy galax- 
ies and PG quasars - derived using three methods and for 
which nuclear 10 /jm fluxes are available (see Table 1 for 
the sample and references for the data). The most direct 
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Table 1 

agns and bh mass methods used for the 10 /zm luminosity vs. bh mass relation. 

Method AGN Rcf 

Reverberation Mapping PG 0026+129*, PG 0052+251*, PG 0804+761*, PG 0844+349*, PG 1211+143* 1 

PG 1226+023*, PG 1229+204*, PG 1307+085*, PG 1351+640, PG 1411+442* 1 

PG 1426+015*, PG 1613+658*, PG 1617+175, PG 1700+51, PG 1704+608 1 

PG 2130+099*, Mkn 335*, Mkn 590*, Mkn 817, NGC 3227*, NGC 4051* 1 

NGC 4151*, NGC 5548*, NGC 7469*, 3C390.3*, Akn 120*, F 9*, IC 4329A*, Mkn 509* 1 

Mkn 279*, NGC 3516*, Mkn 841*, Mkn 766*, NGC 4593* 2 

M BH - cr, relation NGC 1386, NGC 1566, NGC 1667, NGC 2110*, NGC 3185, NGC 3362, NGC 3982 3, 4 

NGC 4388*, NGC 4579*, NGC 5252*, NGC 5273, NGC 5283, NGC 5347 3, 4 

NGC 5695, NGC 5929, NGC 5940, NGC 5953, NGC 6104 , NGC 7682 3, 4 

NGC 7674, Mkn 1, Mkn 3*, Mkn 348*, Mkn 530, Mkn 573, Mkn 744, Mkn 1040 3, 4 

Stellar/Gas kinematics NGC 1386, NGC 4258*, NGC 1068, Circinus, NGC 4945* 5 

NGC 4395* (upper limit) 6 



Note. — The galaxies marked with * were used for the hard X-ray luminosity vs. BH mass relation (Figure 1, right panel). 
References for M BH and method: 1. Kaspi et al. (2000). 2. Compilation of Wandel (2002). 3. Nelson & Whittle (1995). 4. 
Ferrarese & Merritt (2000). 5. Compilation of Moran ct al. (1999). 6. Filippcnko & Ho (2001). 

References for the 10 (im fluxes. PG quasars: Neugebauer et al. (1987) and Sanders ct al. (1989). Sy galaxies: small aperture 
ground-based measurements from Rieke (1978); Maiolino et al. (1995); Krabbe, Boker, & Maiolino (2001); P. Lira (2001, private 
communication), and ISO fluxes from Clavel et al. (2000). 



way is from spatially resolved stellar kinematics (see Ho 
1999 for a review), and measurements of masers in disks 
around massive BHs (see Moran, Grecnhill, & Herrnstein 
1999 for a review). 

In the case of bright active galaxies where the stellar 
emission is swamped by the AGN emission, the best tech- 
nique to derive BH masses is reverberation mapping (see 
a review by Peterson 1993). For sample of AGNs ana- 
lyzed in Kaspi et al. (2000) we have used the BH masses 
(and errors) derived from the mean of the FWHM velocity 
measurements. 

The third and most indirect method uses the empiri- 
cal relation between the BH mass and the stellar velocity 
dispersion (the Mbh - c» relation) found by Ferrarese & 
Merritt (2000): log Mbh = 4.80(±0.54) logcr*-2.9(±1.3). 
The Mbh errors here are dominated by the empirical corre- 
lation uncertainty rather than by the observational errors 
in cr* . In Figure 1 we only plot the errors associated with 
the slope of the relation, that is, A log Mbh = . 54 x log ct* . 

We also explore the 2 — lOkeV hard X-ray luminosity 
vs. BH mass relation for Compton thin AGNs based upon 
the good correlation found between the hard X-ray and the 
10 /im luminosities in Sys and PG quasars (Alonso-Herrero 
et al. 2001, Krabbe et al. 2001). The 2 - lOkeV hard X- 
ray luminosities (corrected for absorption) are from Nan- 
dra & Pounds (1994), Lawson & Turner (1997), Bassani 
et al. (1999), George et al. (2000), Iwasawa et al. (2000), 
and the compilation in Alonso-Herrero et al. (1997). We 
have included a few low luminosity AGNs with available 
hard X-ray fluxes corrected for absorption (although mid- 
infrared fluxes were not available). These are NGC 3031, 
NGC 4594, NGC 6251 (M BH from the compilation of Ho 
2002), NGC 3079 (M BH from Moran et al. 1999) and 
NGC 5194 (M B h derived from the M B h - <?* relation, 
where ct* is from Nelson & Whittle 1995). 

When available we have used distances derived from 
surface brightness fluctuations (Tonry et al. 2001). For 



the rest of the galaxies distances were computed using 
H Q = 75kms" 1 Mpc" 1 . 

3. DISCUSSION AND CONCLUSIONS 

The relations between the BH mass and the 10 /im nu- 
clear luminosity and 2 — 10 keV luminosity for our AGN 
sample are presented in Figure 1 (left and right panel re- 
spectively). The latter relation has been previously an- 
alyzed by Wandel & Mushotzky (1986) and Awaki et al. 
(2001) for a limited number of AGNs. Ho (2002) and Wu & 
Han (2001) have investigated the relationship between the 
radio nuclear and total power and the BH mass and found 
a good correlation for bright AGNs. Ho (2002) however 
finds that weakly active galaxies do not follow the same 
relation. 

In a series of papers, McLeod & Rieke (1995), McLeod, 
Rieke, k Storric-Lombardi (1999) and McLeod & McLeod 
(2001) have demonstrated that there is a relation between 
the nuclear B magnitude of quasars and the near-infrared 
luminosity of the host galaxy (that is, the stellar mass of 
the galaxy). This relation is also present in luminous in- 
frared galaxies when the far-infrared luminosity is used to 
represent the bolometric luminosity (McLeod et al. 1999). 
McLeod & McLeod (2001) interpreted these relations as 
the maximum nuclear luminosity possible in a host galaxy 
of a given near- infrared luminosity (mass). In view of this 
relation and the good correlation between the BH mass 
and the bulge luminosity of the host galaxy it would not 
be unreasonable to expect a relation between the BH mass 
and the AGN luminosity represented by the 10 /zm or the 
hard X-ray luminosities, the latter only for Compton thin 
objects. 

In the hard X-ray vs. Mbh relation the lines (riEdd = 
0.01 and r^Edd = 1) are the Eddington ratios computed 
in Awaki et al. (2001, their figure 3) assuming a bolo- 
metric correction of = 27.2 (note that these authors 
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Fig. 1. — Left panel: Relation between the 10 fim nuclear luminosity and the BH mass. Filled stars are PG quasars, X's are Sy Is, open 
triangles are NLSls, and filled squares are Sy 2s. In the on-line version of this figure the data are color coded according to the method used 
for determining the BH mass: green from stellar/gas kinematics, blue from reverberation mapping and red from the A/bh ~~ CT * relation. Right 
panel: same as left panel but for the hard X-ray (2 — 10 keV) luminosity and BH mass. 



pointed out that such bolometric correction may be lower 
for low luminosity AGNs). All the objects located to the 
right of the rjEdd = 0.01 line are low-luminosity AGNs. Al- 
though this is a good relation, it is only useful for Compton 
thin galaxies (column densities Ah < 10 24 cm~ 2 ). Unfor- 
tunately for a significant fraction of type 2 objects the 
2 — 10 keV X-ray emission is not transmitted through the 
torus (e.g., Bassani et al. 1999), and the 2 — lOkeV lumi- 
nosities cannot be used to trace the obscured BH growth. 

The scatter in the relations may be accounted for by effi- 
ciency variations in the accretion process in different kinds 
of AGN and/or evolutionary effects. The evolution will be 
relevant if the accretion rate is higher in the early growing 
stage of an AGN, as proposed by Mathur, Kuraszkiewicz, 
& Czerny (2001) to explain the properties of narrow line 
Sy 1 galaxies (NLSls). This class of galaxies would repre- 
sent an early stage of the AGN evolution based upon their 
lower BH to bulge mass ratios. The four NLSls (shown as 
open triangles in Figure 1) in our sample show a similar 
behavior to the other Sy galaxies in both panels of Fig- 
ure 1, if only with slightly higher accretion rates than the 
average. However the small number of NLSyls does not 
allow us to reach a firm conclusion. 

Another factor that may introduce some scatter in the 
mid-infrared vs. BH mass relation could be the presence 
of nuclear star formation in Sy 2 galaxies. If star forma- 
tion were present (especially for those galaxies observed 
with the relatively large ISO aperture), we would expect 
an excess of 10 fim emission for a given BH mass. We 
do not find such excesses in our sample of galaxies. On 
the other hand, if the optical depths were high for type 
2 objects even at 10 /xm, then for a given BH mass they 
should lie below type 1 objects. Type 1 objects do not 

5 For reference in this model a BH with an accretion rate of M 
L Edd = 1.5 X 10 38 ^™. (i n erg s" 1 ) L Edd . 



generally lie above type 2 objects in the 10 /Ltm luminosity 
vs. BH mass diagram, suggesting that the reprocessing of 
the 10 /im nuclear emission is not severely affected by ge- 
ometric and optical depth effects (see also Alonso-Herrero 
et al. 2001). 

Hosokawa et al. (2001) and Kawaguchi, Shimura, & Mi- 
neshige (2001) constructed a disk-corona model to account 
for the X-ray through infrared spectral energy distribu- 
tions (SEDs) of quasars and Sy galaxies. These models 
demonstrate that the AGN SEDs are sensitive to both the 
BH mass (Mbh) and accretion rate (M). In the optical, 
UV and soft X-ray spectral ranges decreasing M and/or in- 
creasing A/bh will move the peak of the SED toward longer 
wavelengths (i.e., figure 1 in Hosokawa et al. 2001). On 
the other hand, the hard X-ray and mid-infrared spectral 
shapes remain approximately constant for varying Mbh 
and M. This is expected, as in this model the infrared 
bump is produced by dust heated by the AGN, and mod- 
elled with constant spectral indices, regardless of Mbh or 
M. The infrared power varies so that the luminosity ra- 
tio of the big blue bump to the infrared bump remains 
constant. 

The predicted 10 //m luminosity as a function of Mbh 
from Hosokawa et al. (2001) models is shown in Figure 1 
(left panel) for two accretion rates 5 M = 5 — 0.5 x L^dd/c 2 . 
The third accretion rate M = 0.05 x L^dd/c 2 curve has 
been extrapolated from the above mass accretion rates. 
The 10 /jm luminosities of PG quasars are well reproduced 
with accretion rates of M = 0.5 - 5 x ^Edd/c 2 (as found 
also by Kawaguchi et al. 2001 from the optical to X-ray 
properties), whereas Sy galaxies in general require lower 

= 12 X LEdd/c 2 w1 11 represent a disk with Eddington luminosity, 
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Fig. 2. — Predictions of Hosokawa et al. (2001) disk corona model for the ratio of the B-band, J-band (1.2 fim), JV-band (10 fim) and 
hard X-ray (2 — 10 keV) luminosities to the bolometric luminosity of the AGN as a function of Mbh at 2 = (left panel) and z = 3 in the 
observer's rest-frame (right panel). The accretion rate is M = 5 X -LEdd/c 2 (thin lines, red lines in the color version) and M = 0.5 X i^Edd/c 2 
(thick lines, black lines in the color version). 



accretion rates (see also e.g., Su & Malkan 1989). 

Figure 2 shows the ratio of the optical -B-band, near- 
infrared J-band, 10 /Ltm and hard X-ray to the AGN bolo- 
metric luminosity as a function of A/bh at two different 
redshifts z — (left panel) and z = 3 (in the observer's 
rest-frame, right panel) from Hosokawa et al. (2001) model 
predictions. The accretion rate are M = 5.0 x Z/Edd/c 2 and 
M — 0.5 x Z/Edd/c 2 - From this figure it is clear that the ra- 
tios of the 10 /im and hard X-ray luminosities to the AGN 
bolometric luminosity remain approximately constant for 
increasing AJbh at a given M over redshifts of z = — 3. 
Moreover, the ratio of the 10 /im luminosity to the AGN 
bolometric luminosity is also almost constant for the two 
accretion rates. This shows that the 10 /zm luminosites 
can be used to trace the BH growth by comparing the 
mid-infrared luminosity functions at different epochs. 

The optical to bolometric luminosity ratio varies signifi- 
cantly as a function of the BH mass for redshifts z — and 
z = 3, and the two accretion rates. Although the J-band 
to bolometric luminosity ratio does only vary slightly with 
Mbh at z = and M = 5.0 x Z-Edd/c 2 , it does signifi- 
cantly at z = 3 for both accretion rates. Leaving aside the 
unavoidable obscuration (viewing angle) effects of type 2 
AGNs, the complex dependence of the UV, optical and 
near-infrared SEDs with Mbh and M would make it more 
intricate the interpretation of the time evolution of lumi- 
nosity functions to infer the BH growth. 

There is now growing evidence for a missing population 
of obscured type 2 AGNs, necessary to reproduce the hard 
X-ray background (e.g., Fabian & Iwasawa 1999), as well 
as to reconcile the observed local comoving density of BHs 
with that predicted from the QSO luminosity function at 



z = 3 (Haehnelt, Natarajan, & Rees 1998). This indicates 
that a significant fraction of the accretion by BHs may be 
obscured by dust. The two relations discussed here sug- 
gest that the hard X-ray luminosity and more importantly, 
the mid-infrared luminosity can be effective tools to probe 
the processes taking place in obscured nuclei of galaxies. 
These relations can be used to conduct a detailed census 
of the BH masses over a range of redshifts, and measure 
directly their growth rate as long as an estimate of the ac- 
cretion rate is available. The present X-ray space missions 
are already producing results (e.g., Brandt et al. 2002 
and references therein). The launch of SIRTF will provide 
further capabilities to obtain Mbh for obscured AGNs at 
redshifts of up to z ~ 4. 

Summarizing, for a heterogeneous sample of local AGNs 
we find: (i) the Mbh correlates with the 10 /im nuclear 
luminosity in local Sy galaxies and PG quasars; (ii) the 
available data suggest no significant differences between 
type 1 and type 2 objects, implying that the reprocessing 
of the 10 /jm nuclear emission is not severely affected by 
geometric and optical depth effects; (iii) there is a good 
relation between the AZbh and the 2 — lOkeV hard X- 
ray luminosity, but only for Compton thin galaxies; (iv) 
the Mbh vs. L(10zxm) relation will allow to estimate the 
obscured BH mass growth over cosmological times using 
future observations with SIRTF. 
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the Jet Propulsion Laboratory. 
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